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Reply to the comment by L. Kofman, A. Linde and A. A. Star obinsky (hep-ph/9608341 ) to our 
article "Analytic and Numerical Study of Preheating Dynamics" (hep-ph/960820E 



Preheating is the stage of profuse production of out of 
equilibrium particles either by parametric amplification 
of quantum fluctuations or by the growth of spinodal 
instabilities prior to the reheating stage in inflationary 
cosmologies. Such a mechanism was recognized to be 
more efficient than single particle decay in terms of 
reheating the post inflationary universe. 

The origin of such a mechanism was traced to para- 
metric amplification of quantum fluctuations by the evo- 
lution of the inflaton zero mode. In a previous article 
H we studied such processes, offering a detailed numeri- 
cal analysis in several approximation schemes. We found 
some very interesting new phenomena, especially in the 
case of broken symmetry, where, for "slow-roll" initial 
conditions, the zero mode is driven back to the origin 
and most of the energy has been transferred to the non- 
equilibrium fluctuations. The quantitative features, such 
as time scales for the preheating stage, the rate of parti- 
cle production and the total number of particles produced 
are strongly model dependent. 

The possibility of symmetry restoration via this non- 
equilibrium fluctuations was put forth in reference ||, 
wherein such a mechanism was offered as an explanation 
for the unusual behavior found in ||. In a recent article 
H we provided a more thorough analysis of these new 
phenomena, focusing on the description in Minkowski 
space-time, within the large N approximation. This ap- 
proximation provides a consistent, renormalizable, en- 
ergy conserving scheme that permits an accurate numer- 
ical analysis of the evolution equations. 

In this article we recognized that approximate schemes 
using the Mathieu equation to describe the stage of para- 
metric amplification, though qualitatively correct in the 
sense that a solid intuition on the processes can be ob- 
tained from such an approximation, are quantitatively in 
error. We also recognized that the possibility of symme- 
try restoration is not borne out in the broken symmetry 
phase, ruling out the potential explanations proposed in 
ref. Q for our previous results. In our recent article we 
have provided an analytic and numerical study for this 
phenomenon and explained its physical reasons. 

In a recent comment to our article L. Kofman, A. 
Linde and A. A. Starobinsky (KLS) present a criticism 



of our results. Which we analyze within the context of 
our article below. 



I. THE SETTING: 

As mentioned above, the phenomenon of preheating 
and reheating is strongly model dependent, as well as 
highly non-perturbative in nature. Our study focuses on 
the description of this phenomenon in Minkowski space 
to provide a detailed and deeper understanding within 
a simpler setting. The non-perturbative nature of this 
phenomenon requires that a consistent approximation 
scheme be invoked, and we used the large N expansion 
which is a non-perturbative scheme that can be consis- 
tently improved, unlike the Hartree approximation. 

In this approximation, the dynamics, including backre- 
action effects, can be summarized in the evolution equa- 
tions for the expectation value of the zero mode and the 
fluctuations. We analyzed two different cases: unbroken 
and broken symmetry. The comment by KLS introduces 
expansion the of the universe, which certainly does not 
apply to our study. 



A. Unbroken symmetry: Lame vs. Mathieu 

In the unbroken symmetry case we examined the case 
of very large energy density initially stored in the zero 
mode. In the large N limit the field expectation value 
cj)(t) =< $>(x,t) > (with < • • • > being the expectation 
value in the translational invariant but non-equilibrium 
quantum state) obeys the evolution equation 



4>(t)+M 2 [$] 0(t) = O, 



(1) 



while the mode functions driving the quantum fluctua- 
tions obey: 



Xfc 



(k 2 + M 2 [<S>}) X k(t) = 



(2) 



The effective "mass" term M 2 [<&] includes the evolu- 
tion of the zero mode and the back-reaction effects of 
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the mode themselves through the quantum fluctuations. 
Therefore this is a highly non-linear integro-differential 
problem. We emphasized that the full backreaction prob- 
lem is energy conserving and provided a detailed account 
of the renormalization aspects. Note that none of these 
issues were taken up by KLS in their analysis. 

In the initial stage, all the energy is assumed to be in 
the zero mode of the field. This corresponds to choos- 
ing as initial state the "vacuum" for the k ^ modes. 
For very weak coupling the back-reaction effect will be 
negligible up to a time we call the preheating time; we 
will estimate this time below. Thus for the early stages 
of the dynamics and for very weak coupling one can ap- 
proximate 
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(3) 



Now, using the classical oscillating behaviour of (f>(t) 
one is led by eq.(||) to an effective mass that oscillates 
in time. In this approximation (which is indeed very 
good for weak coupling and during the preheating stage, 
see ||), eq.(0) exhibits parametric resonance as noticed 
first in rcf Jy. Namely, there are allowed and forbidden 
bands in k 2 . The modes within the forbidden bands grow 
exponentially, whereas those in the allowed bands remain 
bounded in absolute value. The growth of the modes in 
the forbidden bands is interpreted as the production of 
particles are with these momenta and the production rate 
is determined by the imaginary part of the Floquet index. 
This exponential, non-perturbative growth of quantum 
fluctuations will lead to strong backreaction effects once 
the contribution from the quantum fluctuations to M 2 [$] 
becomes of the same order as the tree-level terms. When 
this happens the approximation ([?]) breaks down, and the 
energy transferred from the zero mode to the produced 
particles is a large fraction of the initial zero mode energy. 

In approximations that do not include the backreac- 
tion effects there is infinite particle production, since the 
effect of draining energy from the zero mode is not taken 
into account. The full backreaction problem and the ap- 
proximation used in our work maintains energy conserva- 
tion, and thus displays the shutting down of the particle 
production when the back-reaction becomes important; 
this is the end of the preheating stage. 

Now, in order to compute quantitatively the number 
of particles produced, the behavior of the zero mode is 
required. In ref. [jl]]|] <fi(t) is approximated by a cosine 
in the calculations. The mode equations then become 
the Mathieu equation. As shown in ref. ||, the exact 
classical solution is actually a cn Jacobi function. The 
difference between this and a trigonometric function is 
profound. 

Let us now compare the results from the exact mode 
solutions obtained in ref. Q with the Mathieu equation 
approximation to it. In units where m 2 = 1 and setting 

i](t) = <f>(t), one finds 
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where cn stands for the Jacobi cosine and we choose for 
initial conditions 7/(0) = r/o , 7/(0) = 0. 

Inserting this form for 7/(r) in eqs.(E|) and (0) yields 
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Xk(t) = o 



(5) 



This is the Lame equation for a particular value of the 
coefficients that make it solvable in terms of Jacobi func- 
tions (see and references therein). As shown in ref. 
this equation has only one forbidden band for posi- 
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tive k 2 , which runs from k 2 — to k 2 = One can 
choose Floquet solutions of eq. (p|) fulfiling the relation 



U k (t + 2uj) = e lF ^ U k (t), 



(6) 



where the Floquet indices F(k) are independent of t. In 
the forbidden band the F(k) posses an imaginary part. 
Their exact form results || 

F(k) = -2iK{k) Z(2K(k) v)+tt 

where Z(u) is the Jacobi zeta function (see ^ and refer- 
ences therein) and v is a function of fc in the forbidden 
band defined by 



^j=cn(2K(k)v,k) , < v < -. 
v 2 2 



(7) 



All these elliptic functions posses fastly convergent ex- 
pansions in powers of the elliptic nome 



-KK'(k)/K(k) 



(8) 



Since < k < \j\[2 [see eq.(|])], we have 
< q < e" 7r = 0.0432139 .... 
Then, 

F(k) = 4i 7r q sin(27rw) [l + 2 q cos 2ttv + 0(q 2 )] + tt . 

(9) 

The imaginary part of this function has a maximum at 
k = fci = | r/o (1 - q) + 0{q 2 ) where § 



T = ImF{ki) = 4 tt q + 0(q 3 ) 



(10) 



This simple formula gives the maximun of the imaginary 
part of the Floquet index in the forbidden band with a 
precision better than 8. 10~ 5 . q can be expressed in terms 
of 7/0 as follows M 



V = 77 



1 (l + r, 2 )^-(l + V 2 /2)^ 

2 (l + r/ 2 )V4 + (i + r? 2/2)V4 
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with an error smaller than ~ 1CP 7 . 

Let us now proceed to the Mathieu equation analysis 
of this problem. The cn Jacobi function can be expanded 
as (see M and references therein) 

cn(z, k) = (1 — q) cos(l — 4g)z + q cos 3z + 0(q 2 ) . 

To zeroth order in q we have 

1 + cos(2tyT + r7 2 ) 



v(tf = ^ 



0{q) 



and 2ui = n/ yl + rjQ+0(q). Under such approximations 
eq.da) becomes the Mathieu equation M 



dz 



-j^ + (a - 2q cos 2z) y{z) 



where 



a = 1 
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(ii) 



4(77§ + 1) 



and z = yjrfc + 1 t. Notice that < q < 1/4 in the 
present case. Eq.(O) posses an infinite number of for- 
bidden bands for > 0. The lower and upper band 
edges for the first band are given by M 
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These values must be compared with the exact result for 
the Lame equation given by (||) : k 2 n j = , k 2 up = 
Although the width of the band is well approximated by 
the Mathieu equation, its absolute position is not. The 
numerical values of the maximum of the imaginary part 
of the Floquet index are given in Table I and compared 
with the exact values from eq.(lO). 





Lame 


*F Mathieu 


%error 


l 


0.2258 . . . 


0.20... 


13% 


4 


0.4985 . . . 


0.37... 


35% 


770 — > 00 


4ire- w = 0.5430 . . . 


0.39... 


39% 



TABLE I. The maximum of the imaginary part of the 
Floquet index T for the Lame equation and for its Mathieu 
approximation. 



We see that the Mathieu approximation underesti- 
mates the exact result by a fraction ranging from 13 % to 
39%. The second forbidden band in the Mathieu equa- 
tion yields Tuathieu = 0.086 . . . for 770 — > 00. This must 
be compared with Thame = corresponding to the fact 
that there is only one forbidden band in the Lame equa- 
tion. 

In ref. [Q, a large discrepancy between Lame and 
Mathieu Floquet indices has been reported within a dif- 
ferent approximation scheme. 

It is worth mentioning that differences in the Floquet 
indices such as those displayed , enter in the exponent. 
In the large N approximation we see that the discrepancy 
between Mathieu and Lame is very large and therefore 
cause very large errors in quantitative estimates of par- 
ticle production and preheating time. 

Although in some cases as found by KLS the Mathieu 
equation gives a reasonable estimate, in some other cases 
it clearly does not. This was the point raised in our 
article, each case must be treated in its own right. 

For example, the number of particles produced during 
reheating is of the order of the exponential of 2T times 
the reheating time in units of 7r/ yl + ?/q . An error of 
25% in T means an error of 25% in the exponent, so that 
one would find 10 9 instead of 10 12 . 

The Mathieu equation approximation would be exact 
in absence of the A$ 4 inflaton self-coupling. That is, for 
the classical potential 0] 



V = -m 2 <$> 2 



ga 2 ® 2 



(12) 



one can consider as classical solution $(rj) = 
$0 cos(mi), (j = 0. However, the potential ( |l^ ) is unsta- 
ble under renormalization (a $ 4 counterterm is needed 
from the one-loop level). Hence, the A = choice is a 
fine-tuning not protected by any symmetry. 

The mode equations (||) apply to the selfcoupled A $ 4 
scalar field. Models for reheating usually contain at least 
two fields: the inflaton and a lighter field u{x) in which 
the inflaton decays. For a 5 a 2 $ 2 coupling, the mode 
equations for the a field take the form ||,|],^| (adapted 
to Minkowski space which is the focus of our article) 



Vk + [k 2 + 



V k (t) = 



(13) 



A new dimensionless parameter j appears here. Neglect- 
ing the a and $ backreaction, we have 



2 {t) 



(14) 

In ref. we show that abundant particle production 

(app ropriate for rehea ting) occurs even for g = A. 

In hcp-ph/960834l| it is stated that "the main subject 
of investigation of |3| was the theory m 2 4> 2 /2 + A^> 4 /4. 
There are two main regimes there: (f> <^ m/\/X and 
cj) ^> m/y/X. The authors of [|| do not make any distinc- 
tion between these two regimes because they neglect ex- 
pansion of the universe. However, from [|| it follows that 
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in expanding universe there is no parametric resonance at 
all for 4> <C m/yX Therefore all results of || related to 
parametric resonance in this regime do not give a correct 
description of reheating in the theory m 2 <fr 2 /2 + A</> 4 /4." 

Firstly, parametric resonance and universe expansion 
are independent processes that desserve separate inves- 
tigation. Secondly, in our investigations in Minkowski 
spacetime Q we provide analytic expressions which hold 
both for small and large amplitudes. One explicitly sees 
there that the parametric resonance dies in the limit of 
small amplitude. 

Moreover, in ref. Q we studied the small amplitud case 
rj << 1, in Minkowski space and compared explicitly with 
the one-loop and linear relaxation (single particle decay) 
results showing in detail how parametric amplification for 
large amplitudes merges with linear relaxation at small 
amplitudes. 



Therefore, the above statement in hcp-ph/9608341 
cannot be applied to ref. |(| neither to our previous works. 

Eqs.(^3|)-(pT[) become Lame equations when r?(i) is ap- 
proximated by the classical solution in Minkowski space- 
time given by (Q). This Lame equation is solvable in 
closed form when the couplings g and A are related as 
follows 



2.9 ( 
- = n(n 



1) , n=l,2,3,. 



In those cases there are n forbidden bands for k 2 > 0. 
The Lame equation exhibits an infinite number of for- 
bidden bands for generic values of j. The Mathieu and 
WKB approximations have also been applied in the non- 
exactly solvable cases Q]. However, as the above analy- 
sis shows, (see Table I) these results cannot be trusted 
quantitatively. The only available precise method con- 
sists on accurate numerical calculations as those of ref. 
HD (where the precision was at least 10~ 6 ). 

As soon as the quantum fluctuations grow and cease 
to be negligible compared with the the classical contribu- 
tion ([l4]), all the approximations discussed so far (Lame, 
Mathieu, etc.) break down. This time is the so-called 
preheating time t re h H ■ One can estimate t re h by equat- 
ing the zero mode energy ( |l4| ) with the estimation of the 
quantum fluctuations derived from the unstable Floquct 
modes |^J^]. Such estimation yields accurate estimates 
when the Lame Floquet indices are used 0. However, 
as emphasized before, because of the non-perturbative 
nature of this time scale « ln[l/A] Q, differences in the 
Floquct indices stemming from a Mathieu equation ap- 
proximation to a Lame equation lead to severe errors in 
the quantitative estimate of such a time scale. 

Although a fairly accurate estimate of t re h can be ob- 
tained via the Lame equation, with the exact classical 
solution for the zero mode, in order to compute physical 
magnitudes beyond t re h, one must solve self-consistently 
the field equations including the back reaction. Clearly 
this requires a numerical treatment. In ref. |3|,[3| this is 
done for the N — > 00 limit and in ref. |^] to one-loop or- 
der. Such a study led to a very clear physical description 



of the non-equilibrium gas of created particles, and its 
equation of state ||. 

Thus, at the expense of re- iterating our conclusions 
and comments in our article, a Mathieu equation approx- 
imation provides a qualitative, intuitive description of 
parametric amplification, and the non-perturbative phe- 
nomena associated with particle production out of equi- 
librium. However, for a quantitative and trustworthy 
estimate of the physical time scales and production rates 
one must study the proper problem in its full complexity. 
Differences in Floquet indices, propagate exponentially 
making the qualitative estimates based on Mathieu un- 
reliable whenever the Mathieu equation is not the proper 
description. Ultimately one must resort to a full numeri- 
cal scheme to study the dynamics at long times, making 
sure that the approximations involved maintain energy 
conservation (or covariant conservation) . Given that this 
phenomenon is not universal and strongly model depen- 
dent, one equation does not fit all the different situations 
and a careful and consistent analysis of each particular 
scenario is needed. 



B. Broken Symmetry: 

Our study of the dynamics in the broken symme- 
try case in connection with the possibility of symme- 
try restoration as advocated by KLS was sparked by the 
following statement in their article "Nonthermal Phase 
Transitions..." by KLS, Phys. Rev. Lett. 76, 1011, 
(1996) § and again by Kofman Q in his review arti- 
cle. In particular, on page 1012, second column, third 
paragraph of this article: "The mechanism of symme- 
try restoration described above is very general; in par- 
ticular, it explains a surprising behavior of oscillations of 
the scalar field found numerically in the 0(./V)-symmetric 
model of Ref. [6]". This reference is to our previous arti- 
cle ||, and there we studied the same situation that we 
re-analyzed in our latest article: Minkowski space-time, 
large N and broken symmetry phase, with the zero mode 
of the field beginning very close to the origin, i.e. with 
"slow- roll" initial conditions. It is precisely this situation 
that we studied in deeper detail in our latest article and 
concluded after a thorough analytic and numerical study 
that there is no symmetry restoration by the quantum 
fluctuations. Energy is conserved and the sum rule 



- 1 + r){oo) 2 + gS(oo) =0 



(15) 



which is nothing but a statement of massless pions and 
a consequence of Goldstone's theorem, is satisfied. The 
final value 77(00) is obtained dynamically, and depends 
on the initial condition that determines the total energy 
of the system which is conserved by the numerics to one 
part in 10 10 which is our numerical accuracy. 

We have provided exhaustive evidence for this behavior 
both analytically and numerically and showed unequivo- 
cally that there is no possibility of symmetry restora- 
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tion as envisaged by KLS and as quoted by || in the 
situation that we studied. 

We have learned privately that Cooper, Habib, Mot- 
tola and Kluger j| have obtained similar conclusions 
in their thorough and independent study of the broken 
symmetry phase, within the same approximation scheme. 
Their results are complementary to ours, Cooper et. al. 
study the strong coupling regime, but the behavior is 
consistent with the results obtained for the weak cou- 
pling case in our article. 

In their comment, KLS seem to agree with our con- 
clusion that in this particular situation studied in our 
articles, when the expectation value is released near the 
potential hill, symmetry restoration does not occur. We 
have pointed out in H and more recently in || that the 
notion of the effective potential is irrelevant for the dy- 
namics and these effects should be understood as a dy- 
namical change of the effective action. The minima of 
the effective action attained dynamically are very differ- 
ent from those of the effective potential. 

In the situation of 'chaotic initial conditions' but with 
a broken symmetry tree level potential, the issue of sym- 
metry breaking is more subtle. In this case the zero mode 
is initially displaced with a large amplitude and very 
high in the potential hill. The total energy density is 
non-perturbatively large. Classically the zero mode will 
undergo oscillatory behavior between the two classical 
turning points, of very large amplitude and the dynam- 
ics will probe both broken symmetry states. Even at the 
classical level the symmetry is respected by the dynam- 
ics in the sense that the time evolution of the zero mode 
samples equally both vacua. 

This situation is reminiscent of finite temperature in 
which case the energy density is finite and above a critical 
temperature the ensemble averages sample both tree level 
vacua with equal probability thus restoring the symme- 
try. In the dynamical case, the "symmetry restoration" 
is just a consequence of the fact that there is a very large 
energy density in the initial state, much larger than the 
top of the tree level potential, thus under the dynamical 
evolution the system samples both vacua equally. 

Thus the criterion for symmetry restoration when the 
tree level potential allows for broken symmetry states is 
that the energy density in the initial state be larger than 
the top of the tree level potential. That is when the am- 
plitude of the zero mode is such that V(r]o) > V(0). In 
this case the dynamics will be very similar to the unbro- 
ken symmetry case, the amplitude of the zero mode will 
damp out, transferring energy to the quantum fluctu- 
ations via parametric amplification, but asymptotically 
oscillating around zero with a fairly large amplitude. 



pages of text plus twenty figures, we provide a detailed 
and exhaustive analysis of the non-equilibrium issues in 
preheating both in the unbroken and broken symmetry 
phases, in Minkowski space-time, within a well defined 
non-perturbative, consistent, renormalizable and energy 
conserving scheme. We have given all the necessary tech- 
nical details including the error estimates in the numeri- 
cal study and included two lengthy appendices with the 
necessary details for the reader to follow all and every 
step of our analysis. 
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II. CONCLUSIONS 



We here a nswered the points raised by KLS in hep- 
ph/9608341| . We refer to our paper H where in over 50 
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